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Abstract 

We study the influence of Gribov copies on gluon and ghost propagators, evaluated nu- 
merically in pure SU{2) lattice gauge theory in the minimal Landau gauge. Simulations are 
done at four different values of /? (namely (3 = 0, 0.8, 1.6 and 2.7) and for volumes up to 16^ 
(up to 24^ at /? = 1.6). For the gluon propagator, Gribov noise seems to be of the order of 
magnitude of the numerical accuracy, even at very small values of the coupling f5. On the 
contrary, for the ghost propagator, Gribov noise is clearly observable for the three values of 
P in the strong-coupling regime. In particular, data corresponding to the minimal Landau 
gauge are always smaller than those obtained in a generic Landau gauge. This result can be 
qualitatively explained. 

Gauge theories, being invariant under local gauge transformations, are systems with redundant 
dynamical variables, which do not represent true dynamical degrees of freedom. The objects of 
interest are not the gauge fields themselves, but rather the classes (orbits) of gauge-related fields. 
The elimination of such redundant gauge degrees of freedom is essential for understanding and 
extracting physical information from these theories. This is usually done by a method called gauge 
fixing, which is based on the assumption that a gauge-fixing condition can be found which uniquely 
determines a representative gauge field on each orbit. It was pointed out by Gribov [|1| that the 
standard gauge-fixing conditions used for perturbative calculations do not in fact fix the gauge 
fields uniquely: for a non-abelian gauge theory, in the Coulomb or in the Landau gauge, there are 
many gauge equivalent configurations satisfying the Coulomb or Landau transversality condition. 
The existence of these Gribov copies does not affect the results from perturbation theory, but their 
elimination could play a crucial role for non-perturbative features of these theories. 

One of the celebrated advantages of lattice gauge theories is that the lattice provides a regular- 
ization which makes the gauge group compact, so that the Gibbs average of any gauge- invariant 
quantity is well-defined and therefore gauge fixing is, in principle, not required. However, because 
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of asymptotic freedom, the continuum limit is the weak-couphng hmit, and a weak-couphng expan- 
sion requires gauge fixing. Moreover, gauge fixing is used in smearing techniques, and is required 
in order to evaluate quark/gluon matrix elements which can be used to extract non-perturbative 
results from Monte Carlo simulations.^ Thus, one is led to consider gauge- dependent quantities 
on the lattice as well. Unfortunately gauge fixing on the lattice is afflicted by the same problem 
of Gribov copies encountered in the continuum case 0, ^, ^ |^. From a numerical point of view 
the existence of Gribov copies in the lattice Landau gauge is not surprising. In fact, this gauge 
condition [0, § is imposed by finding a gauge transformation {g{x)} which brings the functional 

£u[9] = ^^J2T^'[^-9ix)U,ix)g\x + e,)] (1) 

to a minimum.0 This functional can be seen as the action of a spin glass model 0, namely a 
multi-dimensional system characterized by a strong disorder, and therefore it is plausible that one 
finds a large number of local minima for its energy Su- 

Given the appearance of Gribov copies in numerical studies, we need to understand their in- 
fluence on the evaluation of gauge-dependent quantities. In fact, since usually it is not clear how 
an algorithm selects among different Gribov copies and we do not know which kind of bias they 
introduce on the considered quantities {Gribov noise), numerical results using gauge fixing could 
depend on the gauge-fixing algorithm, making their interpretation conceptually difficult. 

Relatively few studies R 0, |10|, O, O, O, KM have been done in order to analyze the infiuence 



of Gribov copies on lattice quantities. In most cases |jTO|, |TT| it was found that Gribov noise is of the 
same order as the numerical accuracy of the simulations, and that it scales down as a pure statistical 
error. However, in the so-called maximally abelian gauge, Gribov noise seems to be quite large and 
to introduce a clear bias on the number of monopoles [^], and on the value of the abelian string 
tension^ Finally, in the case of U{1) lattice gauge theories, it has been shown numerically 

1^, that the photon propagator is strongly affected by Gribov noise and that only averages taken 
on absolute minima of the minimizing functional — the [/(l)-analogue of Su[g] defined in eq. (||) 
— reproduce the theoretical predictions. This suggests that the same could happen for the gluon 
and the ghost propagators. 

In this work we present an extensive study of Gribov noise for gluon and ghost propagators at 
four different values of (3, i.e. (3 = 0,0.8, 1.6 and 2.7. Although these propagators are non-gauge- 
invariant quantities, the study of their infra-red behavior provides a powerful tool for increasing 
our understanding of QCD and, in particular, for gaining insight into the physics of confinement in 
non-abelian gauge theories. Notice that the first three values of (3 are in the strong-coupling regime, 
where the number of Gribov copies is higher and Gribov noise, if present, is probably larger and 
more easily detectable. 



Gluon propagators in Landau gauge have been the subject of several numerical studies P, ITBI, IT? 



However, in these works Gribov copies were not taken into account. In reference [jTT[ Gribov noise 
has been estimated for the electric and the magnetic screening masses — obtained from the behavior 
of the gluon propagator in the deconfined phase — and found to be of the order of magnitude of 



the numerical accuracy. Finally, only one numerical study |18| is available for the ghost propagator. 



^ For a review see [|| . 

^ This definition applies to SU{2) lattice gauge theory in 4 dimensions. Here {U^{x)} is a given (i.e. fixed) 
thermalized lattice configuration. See Section |2| for notation and definitions. 

The string tension, as usually defined, is a gauge-invariant quantity. However, in the maximally abelian gauge, 
the gauge field is projected onto an abelian lattice gauge field, and the string tension constructed in terms of this 
abelian field is not gauge-invariant. For details see p^ . 
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Also in this case Gribov noise was not directly analyzed. However, since the fluctuations of the 
data appear to be very small, the authors concluded that, if present, Gribov noise should be a small 
effect. 



1 Minimal Landau Gauge 

In order to get rid of the problem of spurious gauge copies, Gribov proposed the use of additional 
gauge conditions. In particular he restricted the physical configuration space to the region 

Q = {A : d- A = 0, M[A] > 0} , (2) 

where Ai [A] = —V ■ D [A] is the Faddeev-Popov operator. This region is clearly included in the 
hyperplane F of transverse configurations {i.e. d ■ A = 0) and is delimited by the so-called first 
Gribov horizon, where the smallest eigenvalue^ of the Faddeev-Popov operator is zero. 

Nowadays we know that the region Q is not free of Gribov copies and that conditions stronger 
than those in (H) should be imposed. In particular the physical configuration space has to be 
identified with the so-called fundamental modular region A, which is defined (in the continuum) as 
the set of absolute minima of the functional 



This corresponds to selecting, on a given orbit, the configuration closest to the origin. Of course one 
has to show that this functional possesses a unique absolute minimum on each orbit. This has been 
proven in . From these works we also know that Q is simply the region of local minima of the 
functional -EAif?], namely it includes A, and is not free of Gribov copies. It has also been shown |^ 



that in the interior of the fundamental modular region A the absolute minima are non-degenerate. 
However, on the boundary of A there are degenerate absolute minima, and they have to be identified 
in order to obtain a region truly free of Gribov copies.^ 

To eliminate Gribov copies on the lattice we can define a fundamental modular region also in 
this case. To this end we can look for the absolute minimum of the minimizing functional Su[g] 
{minimal Landau gauge), defined in (^, instead of looking for a generic relative minimum. This 
definition is the lattice analogue of the one used in the continuum [see eq. (|^)]. In fact, it is clear 
that, by choosing the absolute minimum, we select on each gauge orbit the configuration closest 



to the vacuum V^(a;) = IL . It can also be shown |]22|, that the set of all minima — local 



and absolute ones — of the functional Su[g] is the lattice analogue of the region f2 introduced 
by Gribov, namely it is the set of transverse configurations for which the lattice Faddeev-Popov 
operator A4[U] = —V ■ D[U] is non-negative.| 

Of course, on a finite lattice, the existence of an absolute minimum for Su[g] is manifest. In fact 
the gauge orbit is compact and this functional is bounded. However, from a numerical point of view, 
the search for this absolute minimum is a difficult task. The reason is that eq. (|^) represents the 
energy of a nonlinear sigma-model with variables g{x) taking values on a group [in our case SU{2)] 
and, due to this non-linearity, the search for the absolute minimum is highly non-trivial. Moreover, 

The Faddeev-Popov operator has a trivial null eigenvalue, corresponding to a constant eigenvector. 
^ This result does not contradict Singer's no-go theorem since the hypothesis of continuity of the gauge- fixing 
condition is dropped in this approach. 

^ Here V is the lattice four-divergence operator and D[U] is the lattice gauge-covariant derivative. 
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there are no theoretical predictions that allow us to distinguish a generic relative minimum from the 
absolute one, except for looking at all minima of the minimizing functional and checking that the 
(supposed) absolute minimum effectively corresponds to the smallest value of this functional.^ The 
search is therefore a statistical one: enough gauge-fixed gauge copies Cj — for each thermalized 
configuration C — have to be produced so that the probability to get, with a new copy, a value 
of the functional Su[g] smaller than all the previous ones is negligible. In particular, we need a 
statistical criterion in order to decide when to stop this search. Of course, this criterion must 
balance the necessity of producing enough Gribov copies and the need to keep the computational 
work limited to a reasonable amount of CPU time. We have decided to proceed in the following 
way. Given a thermalized configuration C, we start by producing five gauge-fixed gauge copies Cj. 
The copy characterized by the smallest value for the minimizing functional is chosen as a candidate 
for the absolute minimum Cmin- Then we keep doubling the number of gauge- fixed gauge copies Cj. 
After each doubling we check if, in one of the new copies, the minimizing functional gets a value 
smaller than that of the supposed absolute minimum Cmin- In the case of a negative result for this 
check, the search for the absolute minimum is stopped and Cmin is used as a chosen candidate for 
it. In the opposite case, the new gauge-fixed configuration, characterized by the smallest value for 
the minimizing functional, is chosen as a new candidate for the absolute minimum Cmin, and we 
double again the number of Gribov copies. When 320 copies have been produced the search stops 
anyhow, and the configuration which is the present candidate for Cmin is considered as the real 
absolute minimum. So, for any configuration C, a minimum of 10 and a maximum of 320 Gribov 
copies are produced. 

Finally, we consider two different averages: the average considering only the supposed absolute 
minima (denoted by "am"), which should give us the result in the minimal Landau gauge; and the 
average considering only the first gauge-fixed gauge copy generated for each configuration (denoted 
by "fc"). The latter average is the result that we would obtain if Gribov noise were not considered. 

2 Propagators on the Lattice 

We consider a standard Wilson action for SU (2) lattice gauge theory in 4 dimensions with periodic 
boundary conditions. The symbol U^{x) G SU{2) is used for link variables, g{x) G SU{2) for site 
variables, represents a unit vector in the positive fi direction, and V = is the lattice volume. 
(Here we always consider lattices with equal sizes in the four directions.) The gauge field is defined 

as 

(4) 

We also define A^(a;) = (l/2i) Tr[A^(a;) cr"], where a" is a Pauli matrix. Note that A^(a;) approaches 
{l/2)[A^^""'^'^]1^{x) in the continuum limit, where [y4'-^°"*'^]^(a;) is the conventional vector potential. 

If the configuration {U^ (x)} is a stationary point of the functional Su[g], defined in eq. ([l|), then 
the lattice divergence of ^^(x) is null, namely 

(V ■ AT (x) ^ ^» - - e^) = V X, a . (5) 

This is the lattice formulation of the usual Landau gauge-fixing condition in the continuum. By 
summing equation (^ over the components x^ of x with /i 7^ z/, for fixed z/, and using the periodicity 

^ Recently a method to make an extrapolation in the Gribov noise, namely to extrapolate the data in the limit 
of an infinite number of gauge-fixed gauge copies, has been introduced [ p^ . An extension of the present work 
implementing this extrapolation method will be discussed in a future paper. 



M^)^7; U,ix)-Ulix) 
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of the lattice, it is easy to check Q that if the Landau gauge-fixing condition is satisfied then the 
quantities 

Q.M^Y.T.M^) z/ = i,...,4 (6) 

are constant, i.e. independent of x^. 

With this notation the lattice space-time gluon propagator is given by 



D;l{x~y)^{Al{x)K{y))- 



To go to momentum space we can use formula (3.1a) in reference and obtain 

n 2 



1 



12V 



X] 



) 



1 



J2 A^^ix) cos{2nk-x) 



+ 



J2 ^"^(x) sin (27rA; ■ x) 



) 



(7) 

(8) 

(9) 



(Here /z goes from 1 to 4, and k has components A;^ given hj k^N = 0,1,2, . . . , N — 1, where N 
is the lattice size.) Notice that -D(O) is not given by D{k) at k = 0. The difference is due to the 
Landau gauge condition — the continuum- like condition as in equation (^ — which in momentum 
space reads 

Ep(M^;:W ^2Y.sm{7rk,)A';{k) = 0. (10) 

If 7^ we obtain that only three of the four Lorentz components of A^-lk) — and therefore of 
y4"(x) — are independent: this explains the factor 9 (instead of 12) in equation (|^). 
The zero three-momentum-space gluon propagator is defined as P, [if 



(11) 



where invariance under space-time translations is used in order to improve statistics. Q Here i goes 
from 1 to 3 (the three directions orthogonal to the time direction, i.e. i = 4) and 



(12) 



X 



Notice that, if we denote by 0° the zero-momentum component of the gluon field ^"(x), namely 



y r,. 



(where i = 1, 2, 3, 4), then we can write 



V 



and the relation 



E = 



(13) 



(14) 



(15) 



^ The normalization in (|T^) has been chosen so that the first term in equations (|T6|) and (^^ can be directly 
related to the zero- momentum gluon propagator in eq. (p^). 
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holds for any i and a. Thus we obtain, for the zero three-momentum-space gluon propagator, the 
expression 

^(0, ^) = ^ E ( (C)' ) + ^ E ( i^^'^m [^^"^]nt + d) ) (16) 
= $2 + d(o)(o,c/) . (17) 

Clearly, the contribution of the constants 0" to the propagator is considerably enhanced by the 
volume factor in the first term. The effect of the positive constant $^ in eq. (|17D on masses 
obtained in numerical lattice simulations has been recently pointed out in reference |2^. Here 
we evaluate both the total zero three-momentum gluon propagator, as defined in (|l3), and the 
propag ator D(°)(0,rf) defined in (|16D and ([l^). 

Let us observe that a nonzero value for the constants 0^ is related to the use of periodic bound- 
ary conditions and to the finiteness of the volume. In fact, with free boundary conditions, these 
constants are identically null, even on a finite lattice, after the Landau gauge condition is imposed 
Moreover, in the periodic case, it has been proven 123] that these constants must go to 



zero as the volume increases; more exactly, one obtains 0° ^ 1/N as the lattice size goes to 
infinity. This bound is a consequence only of the positivity of the Faddeev-Popov matrix in the 
lattice Landau gauge. Thus, a finite nonzero value of 0° is a lattice artifact due to finite volume and 
translational invariance. Note that this does not imply that should be zero in the infinite volume. 
However, Zwanziger proved ^ ^ that, in the infinite-volume limit, the zero four-momentum 



gluon propagator -D(O), which can be written as0 

^(0) = (18) 

and is clearly related to is less singular at momentum k = than and that, very likely, 
it does vanish as rapidly as k"^. Thus, the infinite- volume limit of is related to the infra-red 
behavior of the gluon propagator, a problem still not settled. 



Finally, following |18|, ^ , we define the space-time ghost propagator as 

irab — / / A ^-1 ' 



G{x-y)6'^'' ^ {(M-') [U]). (19) 



xy 



Going to momentum space and taking the trace in color space, we obtain 

When the lattice Landau gauge is fixed, the Faddeev-Popov matrix is symmetric and semi-positive- 
definite and we can write 



M[U] = -(V - Dp]) = - {Dp] -V) . (21) 

(For an explicit definition of the matrix A^"y[f/] see eq. (B.18) in reference I^Sj.) In particular, 
this matrix is positive-definite in the subspace orthogonal to constant vectors. Therefore, it can 
be inverted by using a standard conjugate-gradient method (CG), provided that we work in this 
subspace. Moreover, since we use this method in order to evaluate the Faddeev-Popov propagator, 

^ Notice that here fi goes from 1 to 4, while in eq. the index i goes from 1 to 3. 
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we need to impose this restriction in a way that preserves translational invariance for the lattice. 
This is not a problem if the source 'ip'^{x) and the initial guess of the solution for the CG- method 

As a source we have used 



have zero constant mode [IS, 27 



6" 



-27rj k-x 



(22) 



where c has a fixed value (1,2 or 3) and k ^ 0. In this way, the condition J2x '4'"'{x) = is 
automatically imposed, and at the same time we save computer time by evaluating "half" of the 
Fourier transform in eq. (pO]). 



3 Numerical Simulations and Results 

To thermalize the gauge configuration {U^{x)} at a fixed value of the coupling P we use a hybrid 
overrelaxed (HOR) algorithm |2^. This algorithm is able to speed up numerical simulations in the 
case of pure lattice gauge theory: in fact, an overall efficiency gain of up to a factor 10 has been 
observed |^9|] [for the SU{2) case] with respect to the heat-bath algorithm, and its dynamic critical 
exponent has been found |3l| to be equal to 1.0(1) [again in the SU(2) case]. The HOR algorithm 



is also easy to be implemented: m micro canonical (or energy-conserving) update sweeps are done, 
followed by one standard local ergodic update (heat-bath sweep) of the lattice. In order to optimize 
the efficiency of the heat-bath code, we implement two different SU{2) generators (methods 1 and 
2 described in Appendix A of [0, with hcutoff = 2). In our case we did not tune the value of m; 
following [^ we fixed m = A^/2. However, for all the pairs {f3,N), we evaluated the integrated 
autocorrelation tim^^ for the Wilson loops W{1, 1) with / = 1,2,4, . . . , A^/2 , and for the Polyakov 
loop P (denoted respectively by Tint,Wi Tint,p)- For all pairs and observables, except for P at 
P = 2.7, we obtained Tint ^ 1- Note that Tint = 0.5 indicates that two successive configurations 
generated in the Monte Carlo simulation are independent. The situation is different at /3 = 2.7, 
where there is clearly a slow mode associated with the Polyakov loop P, a result in agreement with 
the analysis in reference [^. At any rate, also in this case we obtained Tint,p ^ 20 for all lattice 
sizes. Since, for all pairs {jS, N) and for all quantities, the integrated autocorrelation time is much 
smaller than the number of sweeps between two consecutive configurations used for evaluating the 
propagators (see Table 0), we can conclude that these configurations are essentially statistically 
independent. For the lattice 8^ at /? = 2.7 we also made a comparison between the heat-bath and 
the HOR algorithm. We obtained Tint,p = 225 ± 68 in the first case and Tint,p = 12.5 ± 1.2 in the 
second one (see [^ for details). Since, in our case, the deterministic update of a link is almost five 



times faster than the heat-bath one, and "one sweep" of the lattice for the HOR algorithm means 
m = 4 microcanonical sweeps followed by one heat-bath sweep, we gain a real {i.e. CPU-time) 
factor of order 10, in agreement with [^ . 



For the numerical gauge fixing we use the so-called stochastic overrelaxation algorithm 
which has dynamic critical exponent z equal to 1 (see [§6]). In all our simulations we stop the 
gauge fixing when the condition 

[(V-A)'^(x)]' < 10-^2 (23) 



For a definition of dynamic critical exponent see for example ]30| . 

For a definition see pO[ |. To evaluate the integrated autocorrelation time we use an automatic windowing 
procedure ||3^ with two different window factors (6 and 15). We also employ a method Q based on a comparison 
between the naive statistical error with a jack-knife binning error §. In all cases we checked that these three 
estimates are in agreement. 
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is satisfied. This is equivalent to fixing the minimizing functional up to about one part in 10^^, 
and it is sufficient in order to distinguish two different lattice Gribov copies P, |, ^, H . In the final 
gauge-fixed configuration we also evaluate 



Q 



12 N 



Ql Ql 



(24) 



where Qy = ( Y.x^ Qvi^v) , and Qy{xy) has been defined in eq. (P). This quantity should be 
zero when the configuration is gauge-fixed, and it is a good estimator of the quality of the gauge 
fixing. As in reference I^Bj, we found that the stochastic over relaxation algorithm is very efficient 
in fighting critical slowing-down and in making the quantities Qu{xu) constant. 

In Table |l] we report, for each pair {(3,N), the parameters used for the simulations. Overall, we 
have considered about 2000 configurations (gauge orbits) and 40, 000 Gribov copies. In all our runs 
we have started from a randomly chosen lattice gauge configuration. 

From tests conducted at lattice sizes up to 8^, we noticed a strong violation of rotational invari- 
ance, as is expected for small lattices and non-improved actions. Therefore, we decided to evaluate 
the four-momentum-space gluon and ghost propagators by considering only values of k with three 
of the four components equal to zero, namely k = (0, 0, 0, k^). 

Finally, for the ghost propagator at j3 = 1.6, data have been collected only up to the volume 
16^, due to limitations on the memory of the available computers. 

Computations were performed on several IBM RS-6000/250-340 workstations at New York 
University. 



3.1 Gluon Propagator 

In Table ^ we show the results (averages "am" and "fc") for the four-momentum-space gluon 
propagator D{k) as a function of Nk^. The data corresponding to the minimal Landau gauge 
(absolute minima) are in complete agreement, within statistical errors, with those obtained in a 
generic Landau gauge. This happens even at /3 = 0, where the number of Gribov copies is very 
large (see Table p. Clearly for the zero three- momentum-space gluon propagator the result is the 
same: Gribov noise, if present, seems to be irrelevant compared to the statistical fluctuations. This 
is the case for both propagators, namely D(0,d) defined in eq. (|lT]) and D^'^\0,d) defined in eq. 
m and (0). 

However, a few interesting observations can be made about these data. 

First of all, from Table ^we can see that, in the strong-coupling regime, the infra-red behavior 
of the gluon propagator is completely different from that at /? = 2.7, where it grows very fast as the 
momentum goes to zero. In particular, at /5 = 0, this propagator is clearly decreasing as k^ goes 
to zero. At (3 = 0.8 we see an almost constant propagator. Finally, at /3 = L6 the propagator is 
increasing as k^ goes to zero for the lattice size N = 8, while it is constant or slightly decreasing 
with for = 16 and N = 24. Also the behavior of zero-momentum gluon propagator -D(O) as 
a function of the lattice size A^ is very /3-dependent. In fact, at /3 = 0, it decreases as A^ increases, 
and from our data it is not clear if its value would be zero, or a finite nonzero constant, in the limit 
of A^ going to infinity. On the contrary, at /5 = 2.7, this quantity is increasing with the lattice size. 
Finally, at the two intermediate values of f3 we see a very small volume dependence. 

Related to these observations is the fact that, in the strong-coupling regime, the gluon propagator 
D{0,d) is not always positive (see Figure |I|). This is a manifest violation of reflection positivity 
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Finally, looking at Figure ||, it is clear that the positiveness of -D(0, d) at P = 2.7 is only due to 
the large value of the positive constant in eq. (p!7|), as predicted in reference Fig. 1]. 



3.2 Ghost Propagator 

Results for the ghost propagator are reported in Table ^ In this case it can be noticed that the 
data obtained in the minimal Landau gauge (average "am") are constantly smaller than or equal 
to the corresponding "fc"-data.0 Moreover, when the statistics are good, this difference is larger 
than the numerical uncertainties, especially at small values of Nk^. The only case in which this 
does not happen is f3 = 2.7. However, in this case, even for = 16, almost no Gribov copies were 
produced, and therefore we cannot expect a difference between the two sets of data (see Table |I|). 
At any rate, Gribov noise can clearly be observed for the three strong-coupling values of (3. 

This result can be qualitatively explained. In fact, as we said in Section |I], the smallest non- 
trivial eigenvalue Xmin of the Faddeev- Popov operator A^[f/] goes to zero as the first Gribov horizon 
{i.e. the boundary of the region Q) is approached. So, Xmin can be interpreted as a sort of distance 
between the configuration and dQ. Since the fundamental modular region A is included in Q, the 
absolute minimum, which belongs to A, should be "farther" from the the boundary of Q than a 
generic relative minimum. Thus, the absolute minimum configuration should correspond to a value 
of Xmin larger — in average — than the value obtained in a generic relative minimum.|^ This would 
imply [see eq. (|20|)] a ghost propagator smaller (in average) at the absolute minimum, as observed 
in Table ^ 

Recently, Zwanziger proposed a modification of the S'f/(A^)-Yang-Mills action which effec- 
tively constrains the functional integral to the fundamental modular region A in minimal Landau 
gauge. The lattice partition function of this model can be directly evaluated, by the saddle-point 



method, in the thermodynamic limit. From this theory, it is possible to conclude ||2^ that, in four 



Euclidean dimensions, the lattice ghost propagator has a dipole singularity at zero momentum. In 
order to check the infra-red behavior of the ghost propagator we divided G{k) (data in Table ^ 
average "am" ) by the corresponding free propagator 

Gf^'^ik) = [4 ^ stn\7rk^) = l/p\k) . (25) 

In Figure ^ we plot the results for the four different values of (3 and lattice size A^ = 16. (Similar 
results can be obtained for different lattice sizes). From these data it is clear that G{k) diverges 
faster than l/p'^{k) as k goes to zero. In order to analyze this singularity we have tried a fit of the 
data in Table |^ (average "am" ) using the Ansatz 

G{k) = ap-\k) + bp~\k) + ... , (26) 

as in reference [^. In all cases we have found that the fit with b set to zero is very poor. When 
the dipole singularity is included, the fit improves and we have a ^ b. However, for large lattice 
sizes and small momenta, the fit always overestimates the numerical values of the data. This is in 
agreement with results in reference |T^, and suggests a singularity at zero momentum smaller than 
p~'^{k). As a second attempt we tried a fit|3 to the Ansatz 

G{k) = a[p-^{k)Y . (27) 

This is not the case only for f3 = 1.6, TV = 16 and Nk4 — 1. 
^■^ A numerical study of this eigenvalue and its dependence on Gribov copies has been done in reference pTf and 
will be reported in a future paper |Q . Results seem to confirm the scenario described here. 

In both cases the fit has been done using a mathematica program implementing the singular- value-decomposition 
method (see for example pa]). 
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In this case the fit is in general quite good. Finally, in order to probe the small- momentum limit, 
we evaluate the exponent s in eq. (^) using only data corresponding to the two smallest momenta 
available for each lattice size. From the results, reported in Table ^, it can be noticed that s 
decreases as the lattice size increases, while it increases as (3 increases. The last result is also clear 
in Figure |. 



4 Conclusions 

Our data show absence of Gribov noise for the gluon propagator and a nonzero Gribov noise for 
the ghost propagator. In the latter case, the effect is small but clearly detectable for the values of 
P in the strong-coupling region. The fact that this noise is not observable at /3 = 2.7 seems to us to 
be related only to the small volumes considered here. Of course this hypothesis should be checked 
numerically. This is, at the moment, beyond the limits of our computational resources. We stress 
that this is the first case of evidence of Gribov noise in a "truly" non-abelian simulation. To date, 
in fact, Gribov noise was observed only in simulations for the U{1) case [H, [13], or using maximally 



abelian gauge M, 12, 13 



As for the infra-red behavior of these two propagators, the data for the ghost propagator show a 
pole "between" the zeroth-order perturbative behavior p~'^{k) — valid at large momenta — and the 
p~'^{k) singularity predicted in [Q. However, in this case, the volumes used do not really allow an 
analysis of the small- momentum limit, and this analysis is complicated by the fact that the ghost 
propagator cannot be evaluated at zero momentum. For the gluon propagator, at least in the strong- 
coupling regime, the data are quite interesting, showing a propagator decreasing as the momentum 
goes to zero. This anomalous behavior, predicted in |22|, is still observable at /3 = 1.6, if large 



volumes are considered. Related to this behavior is the violation of reflection positivity observable 
in Figure |I[ Finally, our data show that the behavior of the zero three-momentum-space gluon 
propagator is strongly affected by the zero-momentum modes of the gluon field |^ . We think that 



these results deserve a more accurate numerical analysis which will be presented in a separate paper 
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N 


config. 


th. sweeps 


sweeps 


P 


copies 


fc ^ am 


0.0 


4 


400(304) 


10 


10 


0.84 


4220(1127) 


159 


0.0 


6 


200(200) 


20 


10 


0.875 


4060(3121) 


167 


0.0 


8 


200(200) 


10 


10 


0.915 


6120(6071) 


179 


0.0 


10 


100(100) 


20 


10 


0.945 


3850(3840) 


89 


0.0 


12 


80(80) 


20 


10 


0.945 


1940(1936) 


71 


0.0 


14 


20(20) 


20 


10 


0.95 


590(590) 


19 


0.0 


16 


10(10) 


20 


10 


0.96 


510(506) 


9 


0.8 


8 


200(200) 


1000 


100 


0.875 


6950(6627) 


177 


0.8 


12 


22(22) 


1375 


125 


0.9 


1240(1239) 


20 


0.8 


16 


9(9) 


1100 


100 


0.94 


140(140) 


8 


1.6 


8 


100(100) 


1100 


100 


0.865 


1910(1545) 


85 


1.6 


16 


20(20) 


1650 


150 


0.92 


910(905) 


19 


1.6 


24 


7(7) 


1650 


150 


0.91 


290(290) 


6 


2.7 


8 


400(47) 


2200 


200 


0.8 


4020(494) 


17 


2.7 


12 


158(25) 


2200 


200 


0.85 


1600(208) 


11 


2.7 


16 


100(26) 


2200 


200 


0.875 


1000(157) 


13 



Table 1: The pairs {f3,N) used for the simulations, the number of configurations (in brackets, the 
number of configurations for which Gribov copies have been found), the number of sweeps used for 
thermalization, the number of sweeps between two consecutive configurations used for collecting our 
data, the parameter p used by the stochastic overrelaxation algorithm, the total number of Gribov 
copies produced (in brackets, the total number of different Gribov copies), and the number of times 
the first copy produced was not the absolute minimum copy. 
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P 
























0.0 


4 




0.208(0.004} 


0.218(0.003) 


0.226(0.006) 














0.0 


4 


fc 


0.207(0.004) 


0.221(0.003) 


0.229(0.005) 














0.0 


6 


am 


0.185(0.006) 


0.189(0.004) 


0.214(0.005) 


0.219(0.007) 












0.0 


6 


fc 


0.183(0.006) 


0.198(0.004) 


0.215(0.005) 


0.222(0.007) 












0.0 


8 


am 


0.162(0.005) 


0.174(0.004) 


0.196(0.004) 


0.209(0.005) 


0.209(0.006) 










0.0 


8 


fc 


0.160(0.005) 


0.177(0.004) 


0.186(0.004) 


0.209(0.005) 


0.211(0.006) 










0.0 


10 


am 


0.154(0.007) 


0.168(0.005) 


0.195(0.006) 


0.211(0.007) 


0.202(0.007) 


0.22(0.01) 








0.0 


10 


fc 


0.157(0.007) 


0.172(0.006) 


0.196(0.006) 


0.203(0.007) 


0.214(0.007) 


0.22(0.01) 








0.0 


12 


am 


0.160(0.007) 


0.147(0.005) 


0.186(0.007) 


0.194(0.007) 


0.203(0.008) 


0.205(0.007) 


0.21(0.01) 






0.0 


12 


fc 


0.148(0.007) 


0.169(0.006) 


0.181(0.007) 


0.192(0.006) 


0.214(0.008) 


0.210(0.007) 


0.22(0.01) 






0.0 


14 


am 


0.15(0.01) 


0.17(0.01) 


0.18(0.01) 


0.19(0.01) 


0.23(0.02) 


0.23(0.02) 


0.23(0.02) 


0.24(0.03) 




0.0 


14 


fc 


0.15(0.01) 


0.14(0.01) 


0.18(0.02) 


0.19(0.01) 


0.18(0.01) 


0.24(0.02) 


0.25(0.02) 


0.19(0.02) 




0.0 


16 


am 


0.12(0.01) 


0.15(0.02) 


0.17(0.02) 


0.16(0.01) 


0.19(0.02) 


0.21(0.02) 


0.20(0.02) 


0.21(0.02) 


0.20(0.03) 


0.0 


16 


fc 


0.12(0.01) 


0.16(0.02) 


0.15(0.02) 


0.22(0.02) 


0.20(0.02) 


0.22(0.02) 


0.20(0.02) 


0.23(0.03) 


0.21(0.04) 


0.8 


8 


am 


0.268(0.008) 


0.284(0.007) 


0.281(0.007) 


0.273(0.006) 


0.251(0.009) 










0.8 


8 


fc 


0.281(0.008) 


0.293(0.007) 


0.281(0.007) 


0.268(0.006) 


0.26(0.01) 










0.8 


12 


am 


0.26(0.02) 


0.23(0.02) 


0.25(0.02) 


0.30(0.02) 


0.28(0.02) 


0.24(0.01) 


0.28(0.03) 






0.8 


12 


fc 


0.25(0.02) 


0.27(0.02) 


0.26(0.02) 


0.30(0.02) 


0.26(0.02) 


0.25(0.02) 


0.28(0.02) 






0.8 


16 


am 


0.26(0.04) 


0.25(0.02) 


0.27(0.02) 


0.39(0.05) 


0.33(0.04) 


0.25(0.02) 


0.23(0.03) 


0.26(0.03) 


0.25(0.04) 


0.8 


16 


fc 


0.23(0.02) 


0.31(0.04) 


0.30(0.04) 


0.28(0.04) 


0.26(0.03) 


0.27(0.04) 


0.25(0.01) 


0.25(0.03) 


0.28(0.04) 


1.6 


8 


am 


0.64(0.03) 


0.61(0.02) 


0.47(0.02) 


0.31(0.01) 


0.27(0.01) 










1.6 


8 


fc 


0.67(0.03) 


0.61(0.02) 


0.47(0.01) 


0.32(0.01) 


0.28(0.01) 










1.6 


16 


am 


0.50(0.04) 


0.55(0.03) 


0.68(0.04) 


0.49(0.03) 


0.47(0.03) 


0.33(0.02) 


0.32(0.02) 


0.29(0.01) 


0.31(0.03) 


1.6 


16 


fc 


0.59(0.05) 


0.58(0.03) 


0.57(0.04) 


0.52(0.05) 


0.45(0.03) 


0.39(0.03) 


0.31(0.02) 


0.29(0.02) 


0.29(0.04) 


1.6 


24 


am 


0.57(0.04) 


0.7(0.1) 


0.60(0.05) 


0.63(0.04) 


0.50(0.06) 


0.49(0.05) 


0.48(0.05) 


0.40(0.09) 


0.35(0.05) 


1.6 


24 


fc 


0.55(0.04) 


0.60(0.06) 


0.58(0.09) 


0.58(0.08) 


0.5(0.1) 


0.42(0.05) 


0.50(0.09) 


0.53(0.07) 


0.45(0.08) 


2.7 


8 


am 


30.3(0.5) 


1.9(0.1) 


0.274(0.005) 


0.141(0.002) 


0.114(0.003) 










2.7 


8 


fc 


30.2(0.5) 


2.0(0.1) 


0.276(0.005) 


0.142(0.002) 


0.114(0.003) 










2.7 


12 


am 


59(2) 


4.2(0.3) 


0.64(0.02) 


0.275(0.007) 


0.169(0.004) 


0.132(0.003) 


0.113(0.004) 






2.7 


12 


fc 


58(1) 


4.4(0.3) 


0.64(0.02) 


0.273(0.007) 


0.171(0.004) 


0.131(0.003) 


0.114(0.004) 






2.7 


16 


am 


94(3) 


8.2(0.8) 


1.19(0.05) 


0.47(0.01) 


0.253(0.009) 


0.183(0.006) 


0.144(0.005) 


0.123(0.004) 


0.116(0.005) 


2.7 


16 


fc 


93(3) 


8.2(0.7) 


1.17(0.04) 


0.47(0.01) 


0.259(0.009) 


0.183(0.006) 


0.144(0.005) 


0.123(0.004) 


0.114(0.005) 



Table 2: The four- momentum-space gluon propagator D{k) [see eq. (||) and (||)] as a function of 
Nk^. In all our runs we set k = (0,0,0, /C4). Clearly the same value of Nk^ for different lattice 
sizes N does not correspond to the same lattice momentum p4 [see eq. (|T0|)]. Since we use periodic 
boundary conditions, only data for Nk4^ < N/2 are reported here. Two different types of statistics 
are considered: "am" and "fc". (For a definition see the last paragraph in Section [^.) Error bars (in 
brackets) are one standard deviation. To make the table more readable, data for = 24, (3 = 1.6 
and ^4 > 8 are not reported. 
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2 4 
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Figure 1: The zero three-momentum-space gluon propagator D{0,d) [see eq. ([TT|) ] as a function 
of the time separation d. Since we use periodic boundary conditions, only data for d < N/2 are 
reported here. Data correspond to = 14^ at /3 = (□), and V = 12^ at /? = 0.8 (O). Here we 
consider the "fc" statistics (for a definition see the last paragraph in Section |T].) Error bars are one 
standard deviation. 
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Figure 2: Comparison between the zero three-momentum-space gluon propagators D{0, d) (□) and 
D^^\{)^d) {<>) [see eqs. (0) and (|16|) ] as a function of the time separation d. Since we use periodic 
boundary conditions, only data for d < N/2 are reported here. The data, normahzed to one at the 
zero time-separation, correspond ioV = 16^ and (3 = 2.7. Here we consider the "fc" statistics (for 
a definition see the last paragraph in Section ^) Error bars are one standard deviation. 
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Af 


Stat. 


1 


2 


3 


4 


5 


6 


7 


8 


0.0 


4 


am 


2.29(0.03) 


0.894(0.007) 














0.0 


4 


fc 


2.35(0.03) 


0.916(0.008) 














0.0 


6 


am 


5.21(0.05) 


1.287(0.009) 


0.876(0.005) 












0.0 


6 


fc 


5.54(0.06) 


1.34(0.01) 


0.906(0.007) 












0.0 


8 


am 


9.85(0.08) 


2.19(0.02) 


1.088(0.005) 


0.880(0.005) 










0.0 


8 


fc 


10.19(0.08) 


2.23(0.01) 


1.102(0.004) 


0.889(0.004) 










0.0 


10 


am 


16.0(0.1) 


3.48(0.02) 


1.534(0.005) 


1.004(0.003) 


0.877(0.003) 








0.0 


10 


fc 


16.9(0.2) 


3.58(0.02) 


1.564(0.007) 


1.015(0.004) 


0.888(0.003) 








0.0 


12 


am 


24.1(0.2) 


5.19(0.02) 


2.177(0.008) 


1.285(0.004) 


0.960(0.002) 


0.876(0.002) 






0.0 


12 


fc 


24.6(0.2) 


5.28(0.03) 


2.204(0.009) 


1.296(0.004) 


0.968(0.003) 


0.883(0.003) 






0.0 


14 


am 


33.1(0.3) 


7.28(0.04) 


2.99(0.01) 


1.681(0.007) 


1.157(0.004) 


0.938(0.003) 


0.876(0.003) 




0.0 


14 


fc 


34.8(0.3) 


7.54(0.04) 


3.07(0.01) 


1.706(0.006) 


1.170(0.004) 


0.948(0.003) 


0.888(0.004) 




0.0 


16 


am 


44.2(0.7) 


9.9(0.1) 


4.04(0.04) 


2.20(0.02) 


1.434(0.006) 


1.085(0.004) 


0.926(0.004) 


0.879(0.004) 


0.0 


16 


am 


45.5(0.7) 


9.9(0.1) 


4.02(0.04) 


2.18(0.02) 


1.432(0.009) 


1.085(0.006) 


0.922(0.005) 


0.874(0.005) 


0.8 


8 


am 


8.41(0.07) 


1.778(0.008) 


0.861(0.003) 


0.689(0.002) 










0.8 


8 


fc 


8.94(0.08) 


1.843(0.009) 


0.880(0.003) 


0.701(0.003) 










0.8 


12 


am 


20.8(0.1) 


4.42(0.03) 


1.793(0.009) 


1.031(0.004) 


0.761(0.003) 


0.690(0.002) 






0.8 


12 


fc 


22.1(0.6) 


4.53(0.05) 


1.84(0.01) 


1.053(0.005) 


0.772(0.003) 


0.697(0.003) 






0.8 


16 


am 


39.3(0.6) 


8.47(0.08) 


3.37(0.02) 


1.800(0.009) 


1.160(0.006) 


0.866(0.004) 


0.732(0.003) 


0.693(0.003) 


0.8 


16 


fc 


40.6(0.7) 


8.7(0.1) 


3.43(0.03) 


1.82(0.02) 


1.18(0.01) 


0.875(0.007) 


0.738(0.005) 


0.698(0.004) 


1.6 


8 


am 


6.61(0.06) 


1.335(0.006) 


0.632(0.002) 


0.503(0.001) 










1.6 


8 


fc 


7.2(0.1) 


1.38(0.01) 


0.645(0.004) 


0.512(0.003) 










1.6 


16 


am 


32.2(0.5) 


6.66(0.07) 


2.56(0.02) 


1.329(0.005) 


0.845(0.003) 


0.627(0.002) 


0.528(0.001) 


0.500(0.001) 


1.6 


16 


fc 


32.1(0.3) 


6.66(0.04) 


2.56(0.01) 


1.332(0.006) 


0.847(0.003) 


0.628(0.002) 


0.530(0.001) 


0.5004(0.0009) 


2.7 


8 


am 


3.4(0.1) 


0.635(0.002) 


0.3494(0.0004) 


0.2937(0.0003) 










2.7 


8 


fc 


3.4(0.1) 


0.636(0.002) 


0.3495(0.0005) 


0.2938(0.0003) 










2.7 


12 


am 


7.0(0.2) 


1.308(0.005) 


0.610(0.001) 


0.3939(0.0004) 


0.3116(0.0002) 


0.2893(0.0002) 






2.7 


12 


fc 


7.1(0.3) 


1.309(0.005) 


0.610(0.001) 


0.3939(0.0004) 


0.3117(0.0002) 


0.2894(0.0002) 






2.7 


16 


am 


12.7(0.5) 


2.31(0.01) 


1.016(0.002) 


0.6042(0.0006) 


0.4272(0.0003) 


0.3412(0.0002) 


0.3004(0.0001) 


0.2883(0.0001) 


2.7 


16 


fc 


12.9(0.6) 


2.32(0.01) 


1.017(0.002) 


0.6044(0.0006) 


0.4273(0.0003) 


0.3413(0.0002) 


0.3004(0.0001) 


0.2883(0.0001) 



Table 3: The four-momentum-space ghost propagator G{k) [see eq. (^)] as a function of Nk4. 
In all our runs we set k = (0,0,0, /C4). Clearly the same value of Nk4 for different lattice sizes 
N does not correspond to the same lattice momentum [see eq. (|l3)]- Since we use periodic 
boundary conditions, only data for Nk^ < N/2 are reported here. Two different types of statistics 
are considered: "am" and "fc". (For a definition see the last paragraph in Section ||.) Error bars 
(in brackets) are one standard deviation. 
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Figure 3: Plot of the four-momentum-space ghost propagator G{k) [see eq. (pOl)] divided by the 
corresponding free propagator G-^'^'^^^k) [see eq. (^51)] as a function of p^(A;) [see eq. (j^)]- In all our 
runs we set k = (0, 0, 0, k^). Since we use periodic boundary conditions, only data for k^ < 1/2 are 
reported here. The data, normalized to one at the smallest nonzero momentum k^ = 1, correspond 
to 1/ = 16^ at /3 = (□),/? = 0.8 (O), /? = 1.6 (*), /3 = 2.7 (x). Here we consider the "am" statistics 
(for a definition see the last paragraph in Section |I|.) Error bars are one standard deviation. 
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/3 


N 


s 


0.0 


4 


1.35(0.03) 


0.0 


6 


1.27(0.01) 


0.0 


8 


1.23(0.01) 


0.0 


10 


1.187(0.009) 


0.0 


12 


1.166(0.009) 


0.0 


14 


1.135(0.001) 


0.0 


16 


1.11(0.02) 


0.8 


8 


1.27(0.01) 


0.8 


12 


1.18(0.01) 


0.8 


16 


1.13(0.02) 


1.6 


8 


1.30(0.01) 


1.6 


16 


1.17(0.02) 


2.7 


8 


1.36(0.03) 


2.7 


12 


1.27(0.03) 


2.7 


16 


1.27(0.03) 



Table 4: Estimate of the pole of the ghost propagator G{k) [see eq. (^)] as the momentum goes 
to zero. The values of the exponent s [see eq. (0)] are obtained from the data corresponding to 
the two smallest momenta available for each lattice size (Table |], average "am"). Errors bars (in 
brackets) come from propagation of the statistical error on ghost propagator. 
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